set more off
set graphics off
clear all 	
set maxvar 32000
set matsize 800


/*
****************************
****************************
****************************

REPLICATION DO FILE for:
 
Di Leo, Riccardo, Mariaelisa Epifanio, Thomas J. Scotto, and Vera E. Troeger. Forthcoming. "The Motherhood Penalties: Insights from Women in UK Academia" Community, Work & Family doi 10.1080/13668803.2026.2618679doi: 10.1080/13668803.2026.2618679


NOTE: To preserve the anonymity of the respondents, the dataset presents four discrepancies from the one employed for the analysis in the paper. Results may slightly differ from those in the main paper due to the following amendments:

1 - Departments combined (IT/Math, Arch/Eng, Hard Sci -> STEM)
2 - Respondent birth cohorts combined (40s, 50s -> 1950s and Earlier) 
3 - RAE scores coarsened (6 categories -> 3: Low/Medium/High) 
4 - Childbirth decades combined (pre-1990, 1990s -> 1999 and Earlier) 

****************************
****************************
****************************
*/



*CAUTION!!! Local paths have to be adjusted for do-file to work! 
*global root ""
*cd "$root"

*glo tables ""
*glo figures ""

*ssc install outreg2 , replace 
*ssc install copydesc , replace 
*ssc install reghdfe , replace 
*ssc install require , replace
*ssc install ftools , replace


use "motherhood_penalty_dataset" , clear

mkdir "Tables"

**************************************************************************************************************************************************************************************************************************************************** DESCRIPTIVE STATS ************************************************************* 
********************************************************************************************************************************************************************************


################## ################## ################## 
################## TABLE 1 
################## ################## ################## 

ta birth_decade , gen(birth_decade_)
ta curr_responsibmiss , gen(curr_responsibmiss_)
ta curr_dept , gen(curr_dept_)

est clear 

glo controls_sum "curr_gross curr_gross_av_dummy curr_professor curr_permanent curr_startyear curr_satis_reg curr_resfair_reg curr_salfair_reg curr_commute_reg curr_responsibmiss_* curr_russellgroup curr_raescore_2008 curr_acfem_pertotal2015 curr_av_fem_sal13 curr_dept_* birth_decade_* nkids njobs"

eststo mothers: estpost summarize $controls_sum if curr_mothers== 1
eststo non_mothers: estpost summarize $controls_sum if curr_mothers== 0
eststo diff_mothers: estpost ttest $controls_sum , by(curr_mothers)    


esttab mothers non_mothers diff_mothers ///
using "Tables/Table_1", ///
mtitles("Mean (mothers)" "Mean (non-mothers)" "St. Dev. (mothers)" "St. Dev. (non-mothers)" "Count(mothers)" "Count (non-mothers)" "Diff." ) ///
cells("mean(pattern(1 1 0) fmt(2))  sd(pattern(1 1 0) fmt(2))  count(pattern(1 1 0) fmt(0)) b(pattern(0 0 1) fmt(2)) p(pattern(0 0 1) fmt(3))") ///
varl(nkids "N. Children (0-3)" njobs "N. Occupation Stints (1-7)" ///
curr_gross "Current Gross Salary (1-4)" curr_gross_av_dummy "Current Gross Salary above £62.5k (0-1)" curr_professor "Professor (0-1)" curr_permanent "Permanent Contract (0-1)" curr_satis_reg "Satisfied w/ Job (0-1)" curr_resfair_reg "Fair Allocation of Resources wrt. Men (0-1)" curr_salfair_reg "Fair Allocation of Salaries wrt. Men (0-1)" curr_commute_reg "Commuting Time $\leq$ 45 Mins. (0-1)" ///
curr_responsibmiss_1 "Teaching and Research" curr_responsibmiss_2 "Teaching Only" curr_responsibmiss_3 "Research Only" curr_responsibmiss_4 "Others + No Info." ///
curr_startyear_reg "Employment Start Year" ///
curr_russellgroup "Russell Group HEI (1-0)" curr_acfem_pertotal2015 "\% of Female Academics over N. Employed at HEI (2015/16) - sextiles" curr_av_fem_sal13 "Avg. Salary among Female Employees at HEI (2012/13) - sextiles" curr_raescore_2008 "HEI's Research Assessment Exercise (RAE) Score 2008 - sextiles" ///
curr_dept_1 "Architecture and Engineering" curr_dept_2 "Hard Sciences" curr_dept_3 "IT and Mathematics" curr_dept_4 "Arts and Humanities" curr_dept_5 "Clinical Sciences" curr_dept_6 "Business Studies, Economics and Law" curr_dept_7 "Social Sciences" curr_dept_8 "Other" ///
birth_decade_1 "1940s and before" birth_decade_2 "1950s" birth_decade_3 "1960s"  birth_decade_4 "1970s" birth_decade_5 "1980s and after"  birth_decade_6 "Missing Info" ) ///
refcat(curr_gross "\textbf{{Current Job: Role and Salary}}" curr_satis_reg "\textbf{{Current Job: Satisfaction and Attitudes}}" curr_responsibmiss_1 "\textbf{{Current Job: Main Responsibilities}}" curr_russellgroup "\textbf{{Current Job: Academic Info.}}" curr_dept_1 "\textbf{{Academic Discipline}}" birth_decade_1 "\textbf{{Birth Decade}}" , nolabel) ///
stats(N, label("Observations")) ///
replace tex nonotes addnotes("Salary is aggregated to take values 1 'Under £35k', 2 '£35,001-45,000', 3 '£45,001-55,000', 4 'More than 55,001'. Professor is an indicator variable=1 if the respondent reports her current title to be such (else 0). Permanent is an indicator variable=1 if current contract is `permanent' (else 0). All variables are measured at the time of the interview. The different n. of observations are due to attrition into valid responses.")
est clear	

drop birth_decade_* curr_dept_* curr_responsibmiss_*



################## ################## ################## 
################## TABLE 2
################## ################## ################## 

ta young_newbirthdecade , gen(young_newbirthdecade_)

glo youngchild "young_permanent young_excluded_dummy young_relief young_p1_fwe_ml_d young_childcare_dur_dummy young_supportive_before_table young_lesshelp young_morehelp2 young_nohelp2 young_rolehigh young_rolelow young_responsibhigh young_responsiblow young_noresponsibrole young_research_maternity young_teach_maternity young_admin_maternity young_mentor_maternity young_journal_maternity young_russellgroup young_acfem_pertotal2015 young_av_fem_sal13 young_raescore_2008 young_newbirthdecade_1 young_newbirthdecade_2 young_newbirthdecade_3 young_mageatbirth"
*omitted: young_newbirthdecade_4 

est clear

eststo first_birth: estpost summarize $youngchild if knowbirth==1

esttab first_birth ///
using "Tables/Table_2", ///
cells("count(fmt(0)) mean(fmt(2)) sd(fmt(2)) min(fmt(2)) max(fmt(2))") ///
varl(young_permanent "Permanent Contract at Childbirth (0-1)" ///
young_p1_fwe_ml_d "OMP Generosity in FWE (1-6)" young_childcare_dur_dummy "30+ Weekly Childcare Hours (0-1)" young_supportive_before_table "Head of Dept. Support Before Leave (1-3)"  young_excluded_dummy "Felt Excluded upon Return from Leave (0-1)" young_relief "Any Relief Available upon Return from Leave (0-1)" ///
young_lesshelp "Partner Support: 25\% or Below (0-1)" young_morehelp2 "Partner Support: 75\% or Above (0-1)"  young_nohelp2 "Partner Support: No Info (0-1)" young_rolehigh "Partner Employment: High-Rank Academic Job"  young_rolelow "Partner Employment: Other Academic Job" young_responsibhigh "Partner Employment: High-Rank Non-Academic Job" young_responsiblow "Partner Employment: Other Non-Academic Job" young_noresponsibrole "Partner Employment: No Information Provided" ///
young_research_maternity "Research (0-1)" young_teach_maternity  "Ad-hoc teaching (0-1)" young_admin_maternity  "Admin. duties (0-1)" young_conferences_maternity "Conferences (0-1)" young_mentor_maternity  "Mentoring (0-1)" young_journal_maternity "Journal duties (0-1)" ///
young_russellgroup "Russell Group HEI (0-1)" young_acfem_pertotal2015 "\% of Female Academics over N. Employed at HEI (2015/16) - sextiles" young_av_fem_sal13 "Avg. Salary among Female Employees at HEI (2012/13) - sextiles" young_raescore_2008 "HEI's Research Assessment Exercise (RAE) Score 2008 - sextiles" ///
young_mageatbirth "Mother's Age at Childbirth (Anonym.)" young_newbirthdecade_1 "Before 1990" young_newbirthdecade_2 "1990-1999" young_newbirthdecade_3 "2000-2009" young_newbirthdecade_4 "2010-2019" ///
) ///	
refcat(young_lesshelp "\textbf{{Partner Support in the First 2 Years}}" young_rolehigh "\textbf{{Partner Employment at Childbirth}}" young_research_maternity "\textbf{{Activities while on Leave}}" young_russellgroup "\textbf{{Job at Childbirth: Academic Info.}}" young_newbirthdecade_1 "\textbf{{Childbirth Decade}}" , nolabel) ///
stats(N, label("Observations")) ///
replace tex nonotes nonotes
est clear	

drop young_newbirthdecade_*



################## TABLE 3

est clear
eststo full: estpost tab young_p1_fwe_ml_d
la var young_p1_fwe_ml_d "Full Weeks Equivalent (FWE)"

esttab full ///
using "Tables/Table_3", ///
cells("b(pattern(1 1) fmt(0)) pct(pattern(1 1) fmt(2))") ///	
collabels ("Number" "Percent") /// 
replace tex nonotes plain label 





**************************************************************************************************************************************************************************************************************************************************** MOTHERHOOD PENALTY ************************************************************* 
********************************************************************************************************************************************************************************

################## ################## ################## 
################## MAIN 
################## ################## ################## 

################## TABLE 4

est clear

global controls "curr_professor curr_permanent curr_startyear njobs curr_commute curr_russellgroup curr_raescore_2008 curr_acfem_pertotal2015 curr_av_fem_sal13"
global controls2 "curr_startyear njobs curr_commute curr_russellgroup curr_raescore_2008 curr_acfem_pertotal2015 curr_av_fem_sal13"

capture noisily eststo: reghdfe curr_gross $controls curr_mothers , vce(robust) absorb(birth_decade curr_responsibmiss curr_dept)
		estadd loc controls "Yes", replace							
		estadd loc HEI "Yes", replace
		estadd loc dept "Yes", replace							
		
capture noisily eststo: reghdfe curr_professor $controls2 curr_mothers , vce(robust) absorb(birth_decade curr_responsibmiss curr_dept)
		estadd loc controls "Yes", replace							
		estadd loc HEI "Yes", replace
		estadd loc dept "Yes", replace							

esttab est* using "Tables\Table_4", star(* .1 ** .05 *** .01) se b(3) r2 ///
keep(curr_mothers _cons) ///
order(curr_mothers _cons) ///
varl(curr_mothers "Mothers" _cons "Constant") ///
s(controls HEI dept N r2, ///
label("Indiv. Controls" "HEI Controls" "Research Area FEs" "Observations" "R-squared") ///
fmt(0 0 0 0 3)) ///
replace tex nonotes label nobaselevels 



################## TABLE 5

est clear

global controls "curr_professor curr_permanent curr_startyear njobs curr_commute curr_russellgroup curr_raescore_2008 curr_acfem_pertotal2015 curr_av_fem_sal13"
global controls2 "curr_professor curr_permanent curr_startyear njobs curr_russellgroup curr_raescore_2008 curr_acfem_pertotal2015 curr_av_fem_sal13"

capture noisily eststo: reghdfe curr_satis_reg $controls curr_mothers , vce(robust) absorb(birth_decade curr_responsibmiss curr_dept)
		estadd loc controls "Yes", replace							
		estadd loc HEI "Yes", replace
		estadd loc dept "Yes", replace							

capture noisily eststo: reghdfe curr_salfair_reg $controls curr_mothers , vce(robust) absorb(birth_decade curr_responsibmiss curr_dept)
		estadd loc controls "Yes", replace							
		estadd loc HEI "Yes", replace
		estadd loc dept "Yes", replace	
		
capture noisily eststo: reghdfe curr_resfair_reg $controls curr_mothers , vce(robust) absorb(birth_decade curr_responsibmiss curr_dept)
		estadd loc controls "Yes", replace							
		estadd loc HEI "Yes", replace
		estadd loc dept "Yes", replace	
		
esttab est* using "Tables\Table_5", star(* .1 ** .05 *** .01) se b(3) r2 ///
keep(curr_mothers _cons) ///
order(curr_mothers _cons) ///
varl(curr_mothers "Mothers" _cons "Constant") ///
s(controls HEI dept N r2, ///
label("Indiv. Controls" "HEI Controls" "Research Area FEs" "Observations" "R-squared") ///
fmt(0 0 0 0 3)) ///
replace tex nonotes label nobaselevels




################## ################## ################## 
################## APPENDIX C ################## 
################## ################## ################## 

################## TABLE C.1

est clear

global controls "curr_professor curr_permanent curr_startyear njobs curr_commute curr_russellgroup curr_raescore_2008 curr_acfem_pertotal2015 curr_av_fem_sal13"
global controls2 "curr_startyear njobs curr_commute curr_russellgroup curr_raescore_2008 curr_acfem_pertotal2015 curr_av_fem_sal13"

capture noisily eststo: reghdfe curr_gross $controls i.nkids , vce(robust) absorb(birth_decade curr_responsibmiss curr_dept)
		estadd loc controls "Yes", replace							
		estadd loc HEI "Yes", replace
		estadd loc dept "Yes", replace							
		
capture noisily eststo: reghdfe curr_professor $controls2 i.nkids , vce(robust) absorb(birth_decade curr_responsibmiss curr_dept)
		estadd loc controls "Yes", replace							
		estadd loc HEI "Yes", replace
		estadd loc dept "Yes", replace										

esttab est* using "Tables\Table_C1", star(* .1 ** .05 *** .01) se b(3) r2 ///
keep(1.nkids 2.nkids 3.nkids _cons) ///
order(1.nkids 2.nkids 3.nkids _cons) ///
varl(1.nkids "1 Child" 2.nkids "2 Children" 3.nkids "3 Children +" _cons "Constant") ///
s(controls HEI dept N r2, ///
label("Indiv. Controls" "HEI Controls" "Research Area FEs" "Observations" "R-squared") ///
fmt(0 0 0 0 3)) ///
replace tex nonotes label nobaselevels 


################## TABLE C.2

est clear

global controls "curr_professor curr_permanent curr_startyear njobs curr_commute curr_russellgroup curr_raescore_2008 curr_acfem_pertotal2015 curr_av_fem_sal13"
global controls2 "curr_professor curr_permanent curr_startyear njobs curr_russellgroup curr_raescore_2008 curr_acfem_pertotal2015 curr_av_fem_sal13"

capture noisily eststo: reghdfe curr_satis_reg $controls i.nkids , vce(robust) absorb(birth_decade curr_responsibmiss curr_dept)
		estadd loc controls "Yes", replace							
		estadd loc HEI "Yes", replace
		estadd loc dept "Yes", replace							

capture noisily eststo: reghdfe curr_salfair_reg $controls i.nkids , vce(robust) absorb(birth_decade curr_responsibmiss curr_dept)
		estadd loc controls "Yes", replace							
		estadd loc HEI "Yes", replace
		estadd loc dept "Yes", replace	
		
capture noisily eststo: reghdfe curr_resfair_reg $controls i.nkids , vce(robust) absorb(birth_decade curr_responsibmiss curr_dept)
		estadd loc controls "Yes", replace							
		estadd loc HEI "Yes", replace
		estadd loc dept "Yes", replace	
		
esttab est* using "Tables\Table_C2", star(* .1 ** .05 *** .01) se b(3) r2 ///
keep(1.nkids 2.nkids 3.nkids _cons) ///
order(1.nkids 2.nkids 3.nkids _cons) ///
varl(1.nkids "1 Child" 2.nkids "2 Children" 3.nkids "3 Children +" _cons "Constant") ///
s(controls HEI dept N r2, ///
label("Indiv. Controls" "HEI Controls" "Research Area FEs" "Observations" "R-squared") ///
fmt(0 0 0 0 3)) ///
replace tex nonotes label nobaselevels


################## TABLE C.3

est clear

global controls "curr_professor curr_permanent curr_startyear njobs curr_commute curr_russellgroup curr_raescore_2008 curr_acfem_pertotal2015 curr_av_fem_sal13"
global controls2 "curr_startyear njobs curr_commute curr_russellgroup curr_raescore_2008 curr_acfem_pertotal2015 curr_av_fem_sal13"

capture noisily eststo: reghdfe curr_gross $controls nkids , vce(robust) absorb(birth_decade curr_responsibmiss curr_dept)
		estadd loc controls "Yes", replace							
		estadd loc HEI "Yes", replace
		estadd loc dept "Yes", replace							
		
capture noisily eststo: reghdfe curr_professor $controls2 nkids , vce(robust) absorb(birth_decade curr_responsibmiss curr_dept)
		estadd loc controls "Yes", replace							
		estadd loc HEI "Yes", replace
		estadd loc dept "Yes", replace										

esttab est* using "Tables\Table_C3", star(* .1 ** .05 *** .01) se b(3) r2 ///
keep(nkids _cons) ///
order(nkids _cons) ///
varl(nkids "N. Children (0-3)" _cons "Constant") ///
s(controls HEI dept N r2, ///
label("Indiv. Controls" "HEI Controls" "Research Area FEs" "Observations" "R-squared") ///
fmt(0 0 0 0 3)) ///
replace tex nonotes label nobaselevels 


################## TABLE C.4

est clear

global controls "curr_professor curr_permanent curr_startyear njobs curr_commute curr_russellgroup curr_raescore_2008 curr_acfem_pertotal2015 curr_av_fem_sal13"
global controls2 "curr_professor curr_permanent curr_startyear njobs curr_russellgroup curr_raescore_2008 curr_acfem_pertotal2015 curr_av_fem_sal13"

capture noisily eststo: reghdfe curr_satis_reg $controls nkids , vce(robust) absorb(birth_decade curr_responsibmiss curr_dept)
		estadd loc controls "Yes", replace							
		estadd loc HEI "Yes", replace
		estadd loc dept "Yes", replace							

capture noisily eststo: reghdfe curr_salfair_reg $controls nkids , vce(robust) absorb(birth_decade curr_responsibmiss curr_dept)
		estadd loc controls "Yes", replace							
		estadd loc HEI "Yes", replace
		estadd loc dept "Yes", replace	
		
capture noisily eststo: reghdfe curr_resfair_reg $controls nkids , vce(robust) absorb(birth_decade curr_responsibmiss curr_dept)
		estadd loc controls "Yes", replace							
		estadd loc HEI "Yes", replace
		estadd loc dept "Yes", replace	
		
esttab est* using "Tables\Table_C4", star(* .1 ** .05 *** .01) se b(3) r2 ///
keep(nkids _cons) ///
order(nkids _cons) ///
varl(nkids "N. Children (0-3)" _cons "Constant") ///
s(controls HEI dept N r2 , ///
label("Indiv. Controls" "HEI Controls" "Research Area FEs" "Observations" "R-squared") ///
fmt(0 0 0 0 3)) ///
replace tex nonotes label nobaselevels


################## TABLE C.5

est clear

global controls "curr_professor curr_permanent curr_startyear njobs curr_commute curr_russellgroup curr_raescore_2008 curr_acfem_pertotal2015 curr_av_fem_sal13"
global controls2 "curr_startyear njobs curr_commute curr_russellgroup curr_raescore_2008 curr_acfem_pertotal2015 curr_av_fem_sal13"

capture noisily eststo: reghdfe curr_gross $controls multi_nkids , vce(robust) absorb(birth_decade curr_responsibmiss curr_dept)
		estadd loc controls "Yes", replace							
		estadd loc HEI "Yes", replace
		estadd loc dept "Yes", replace							
		
capture noisily eststo: reghdfe curr_professor $controls2 multi_nkids , vce(robust) absorb(birth_decade curr_responsibmiss curr_dept)
		estadd loc controls "Yes", replace							
		estadd loc HEI "Yes", replace
		estadd loc dept "Yes", replace											

esttab est* using "Tables\Table_C5", star(* .1 ** .05 *** .01) se b(3) r2 ///
keep(multi_nkids _cons) ///
order(multi_nkids _cons) ///
varl(multi_nkids "More than 1 Child (One Child)" _cons "Constant") ///
s(controls HEI dept N r2, ///
label("Indiv. Controls" "HEI Controls" "Research Area FEs" "Observations" "R-squared") ///
fmt(0 0 0 0 3)) ///
replace tex nonotes label nobaselevels 


################## TABLE C.6

est clear

global controls "curr_professor curr_permanent curr_startyear njobs curr_commute curr_russellgroup curr_raescore_2008 curr_acfem_pertotal2015 curr_av_fem_sal13"
global controls2 "curr_professor curr_permanent curr_startyear njobs curr_russellgroup curr_raescore_2008 curr_acfem_pertotal2015 curr_av_fem_sal13"

capture noisily eststo: reghdfe curr_satis_reg $controls multi_nkids , vce(robust) absorb(birth_decade curr_responsibmiss curr_dept)
		estadd loc controls "Yes", replace							
		estadd loc HEI "Yes", replace
		estadd loc dept "Yes", replace							

capture noisily eststo: reghdfe curr_salfair_reg $controls multi_nkids , vce(robust) absorb(birth_decade curr_responsibmiss curr_dept)
		estadd loc controls "Yes", replace							
		estadd loc HEI "Yes", replace
		estadd loc dept "Yes", replace	
		
capture noisily eststo: reghdfe curr_resfair_reg $controls multi_nkids , vce(robust) absorb(birth_decade curr_responsibmiss curr_dept)
		estadd loc controls "Yes", replace							
		estadd loc HEI "Yes", replace
		estadd loc dept "Yes", replace	

esttab est* using "Tables\Table_C6", star(* .1 ** .05 *** .01) se b(3) r2 ///
keep(multi_nkids _cons) ///
order(multi_nkids _cons) ///
varl(multi_nkids "More than 1 Child (One Child)" _cons "Constant") ///
s(controls HEI dept N r2, ///
label("Indiv. Controls" "HEI Controls" "Research Area FEs" "Observations" "R-squared") ///
fmt(0 0 0 0 3)) ///
replace tex nonotes label nobaselevels


################## TABLE C.7

est clear

global controls "curr_professor curr_permanent curr_startyear njobs curr_commute curr_russellgroup curr_raescore_2008 curr_acfem_pertotal2015 curr_av_fem_sal13"
global controls2 "curr_startyear njobs curr_commute curr_russellgroup curr_raescore_2008 curr_acfem_pertotal2015 curr_av_fem_sal13"					
		
capture noisily eststo: logit curr_professor $controls2 curr_mothers i.birth_decade i.curr_responsibmiss i.curr_dept , vce(robust)
		estadd loc controls "Yes", replace							
		estadd loc HEI "Yes", replace
		estadd loc dept "Yes", replace							

esttab est* using "Tables\Table_C7", star(* .1 ** .05 *** .01) se b(3) r2 ///
keep(curr_mothers _cons) ///
order(curr_mothers _cons) ///
varl(curr_mothers "Mothers" _cons "Constant") ///
s(controls HEI dept N, ///
label("Indiv. Controls" "HEI Controls" "Research Area FEs" "Observations" "No. of HEIs") ///
fmt(0 0 0 0 0)) ///
replace tex nonotes label nobaselevels 


################## TABLE C.8

est clear

global controls "curr_professor curr_permanent curr_startyear njobs curr_commute curr_russellgroup curr_raescore_2008 curr_acfem_pertotal2015 curr_av_fem_sal13"
global controls2 "curr_professor curr_permanent curr_startyear njobs curr_russellgroup curr_raescore_2008 curr_acfem_pertotal2015 curr_av_fem_sal13"

capture noisily eststo: logit curr_satis_reg $controls curr_mothers i.birth_decade i.curr_responsibmiss i.curr_dept , vce(robust)
		estadd loc controls "Yes", replace							
		estadd loc HEI "Yes", replace
		estadd loc dept "Yes", replace							

capture noisily eststo: logit curr_salfair_reg $controls curr_mothers i.birth_decade i.curr_responsibmiss i.curr_dept , vce(robust)
		estadd loc controls "Yes", replace							
		estadd loc HEI "Yes", replace
		estadd loc dept "Yes", replace	
		
capture noisily eststo: logit curr_resfair_reg $controls curr_mothers i.birth_decade i.curr_responsibmiss i.curr_dept , vce(robust)
		estadd loc controls "Yes", replace							
		estadd loc HEI "Yes", replace
		estadd loc dept "Yes", replace	
		
esttab est* using "Tables\Table_C8", star(* .1 ** .05 *** .01) se b(3) r2 ///
keep(curr_mothers _cons) ///
order(curr_mothers _cons) ///
varl(curr_mothers "Mothers" _cons "Constant") ///
s(controls HEI dept N, ///
label("Indiv. Controls" "HEI Controls" "Research Area FEs" "Observations" "No. of HEIs") ///
fmt(0 0 0 0 0)) ///
replace tex nonotes label nobaselevels


################## TABLE C.9

est clear

global controls "curr_professor curr_permanent curr_startyear njobs curr_commute curr_russellgroup curr_raescore_2008 curr_acfem_pertotal2015 curr_av_fem_sal13"
global controls2 "curr_professor curr_permanent curr_startyear njobs curr_russellgroup curr_raescore_2008 curr_acfem_pertotal2015 curr_av_fem_sal13"

capture noisily eststo: reghdfe curr_satis_reg_alt $controls curr_mothers , vce(robust) absorb(birth_decade curr_responsibmiss curr_dept)
		estadd loc controls "Yes", replace							
		estadd loc HEI "Yes", replace
		estadd loc dept "Yes", replace							

esttab est* using "Tables\Table_C9", star(* .1 ** .05 *** .01) se b(3) r2 ///
keep(curr_mothers _cons) ///
order(curr_mothers _cons) ///
varl(curr_mothers "Mothers" _cons "Constant") ///
s(controls HEI dept N r2, ///
label("Indiv. Controls" "HEI Controls" "Research Area FEs" "Observations" "R-squared") ///
fmt(0 0 0 0 3)) ///
replace tex nonotes label nobaselevels

     

	 
	 

	 
	 
	 
	 
	 
	 
	 
**************************************************************************************************************************************************************************************************************************************************** SURVIVAL MODEL ************************************************************* 
********************************************************************************************************************************************************************************
	 
################## ################## ##################  
################## MAIN
################## ################## ################## 

################## TABLE 6 

est clear

global controls3 "njobs curr_commute curr_startyear curr_russellgroup curr_raescore_2008 curr_acfem_pertotal2015 curr_av_fem_sal13"
global controls3_2 "mothersage_threshold nkids_threshold njobs curr_commute curr_startyear curr_russellgroup curr_raescore_2008 curr_acfem_pertotal2015 curr_av_fem_sal13"

capture noisily eststo: reghdfe curr_professor $controls3 treat_alt , vce(robust) absorb(birth_decade curr_responsibmiss curr_dept)
		estadd loc controls "Yes", replace							
		estadd loc HEI "Yes", replace
		estadd loc dept "Yes", replace							
		estadd loc age2012 "No", replace							
		
capture noisily eststo: reghdfe curr_professor $controls3_2 treat_2nd_alt , vce(robust) absorb(birth_decade curr_responsibmiss curr_dept)
		estadd loc controls "Yes", replace							
		estadd loc HEI "Yes", replace
		estadd loc dept "Yes", replace							
		estadd loc dept "Yes", replace							
		estadd loc age2012 "Yes", replace							

capture noisily eststo: reghdfe curr_gross_av_dummy $controls3 treat_alt , vce(robust) absorb(birth_decade curr_responsibmiss curr_dept)
		estadd loc controls "Yes", replace							
		estadd loc HEI "Yes", replace
		estadd loc dept "Yes", replace							
		estadd loc age2012 "No", replace							
		
capture noisily eststo: reghdfe curr_gross_av_dummy $controls3_2 treat_2nd_alt , vce(robust) absorb(birth_decade curr_responsibmiss curr_dept)
		estadd loc controls "Yes", replace							
		estadd loc HEI "Yes", replace
		estadd loc dept "Yes", replace			
		estadd loc age2012 "Yes", replace							
		
esttab est* using "Tables\Table_6" , star(* .1 ** .05 *** .01) se b(3) r2 ///
keep(treat_alt treat_2nd_alt _cons) ///
order(treat_alt treat_2nd_alt _cons) ///
varl(treat_alt "Becoming Mother (Vs. Not)" treat_2nd_alt "Youngest Child Post-2012 (Vs. Pre-2010)" _cons "Constant") ///
s(controls HEI dept age2012 N r2, ///
label("Indiv. Controls" "HEI Controls" "Research Area FEs" "Age and N. Children by 2012" "Observations" "R-squared") ///
fmt(0 0 0 0 0 3)) ///
replace tex nonotes label nobaselevels 



################## TABLE D.1

est clear

global controls3 "njobs curr_commute curr_startyear curr_russellgroup curr_raescore_2008 curr_acfem_pertotal2015 curr_av_fem_sal13"
global controls3_2 "mothersage_threshold nkids_threshold njobs curr_commute curr_startyear curr_russellgroup curr_raescore_2008 curr_acfem_pertotal2015 curr_av_fem_sal13"

capture noisily eststo: reghdfe curr_gross $controls3 treat_alt , vce(robust) absorb(birth_decade curr_responsibmiss curr_dept)
		estadd loc controls "Yes", replace							
		estadd loc HEI "Yes", replace
		estadd loc dept "Yes", replace							
		estadd loc age2012 "No", replace							
		
capture noisily eststo: reghdfe curr_gross $controls3_2 treat_2nd_alt , vce(robust) absorb(birth_decade curr_responsibmiss curr_dept)
		estadd loc controls "Yes", replace							
		estadd loc HEI "Yes", replace
		estadd loc dept "Yes", replace			
		estadd loc age2012 "Yes", replace							
		
esttab est* using "Tables\Table_D1", star(* .1 ** .05 *** .01) se b(3) r2 ///
keep(treat_alt treat_2nd_alt _cons) ///
order(treat_alt treat_2nd_alt _cons) ///
varl(treat_alt "Becoming Mother (Vs. Not)" treat_2nd_alt "Youngest Child Post-2012 (Vs. Pre-2010)" _cons "Constant") ///
s(controls HEI dept age2012 N r2, ///
label("Indiv. Controls" "HEI Controls" "Research Area FEs" "Age and N. Children by 2012" "Observations" "R-squared") ///
fmt(0 0 0 0 0 3)) ///
replace tex nonotes label nobaselevels 



################## TABLE D.2

est clear

global controls3 "njobs curr_commute curr_startyear curr_russellgroup curr_raescore_2008 curr_acfem_pertotal2015 curr_av_fem_sal13"
global controls3_2 "mothersage_threshold nkids_threshold njobs curr_commute curr_startyear curr_russellgroup curr_raescore_2008 curr_acfem_pertotal2015 curr_av_fem_sal13"

capture noisily eststo: logit curr_professor $controls3 treat_alt i.birth_decade i.curr_responsibmiss i.curr_dept , vce(robust)
		estadd loc controls "Yes", replace							
		estadd loc HEI "Yes", replace
		estadd loc dept "Yes", replace							
		estadd loc age2012 "No", replace							
		
capture noisily eststo: logit curr_professor $controls3_2 treat_2nd_alt i.birth_decade i.curr_responsibmiss i.curr_dept , vce(robust)
		estadd loc controls "Yes", replace							
		estadd loc HEI "Yes", replace
		estadd loc dept "Yes", replace							
		estadd loc dept "Yes", replace							
		estadd loc age2012 "Yes", replace							

capture noisily eststo: logit curr_gross_av_dummy $controls3 treat_alt i.birth_decade i.curr_responsibmiss i.curr_dept , vce(robust)
		estadd loc controls "Yes", replace							
		estadd loc HEI "Yes", replace
		estadd loc dept "Yes", replace							
		estadd loc age2012 "No", replace							
		
capture noisily eststo: logit curr_gross_av_dummy $controls3_2 treat_2nd_alt i.birth_decade i.curr_responsibmiss i.curr_dept , vce(robust)
		estadd loc controls "Yes", replace							
		estadd loc HEI "Yes", replace
		estadd loc dept "Yes", replace			
		estadd loc age2012 "Yes", replace							
		
esttab est* using "Tables\Table_D2", star(* .1 ** .05 *** .01) se b(3) r2 ///
keep(treat_alt treat_2nd_alt _cons) ///
order(treat_alt treat_2nd_alt _cons) ///
varl(treat_alt "Becoming Mother (Vs. Not)" treat_2nd_alt "Youngest Child Post-2012 (Vs. Pre-2010)" _cons "Constant") ///
s(controls HEI dept age2012 N , ///
label("Indiv. Controls" "HEI Controls" "Research Area FEs" "Age and N. Children by 2012" "Observations") ///
fmt(0 0 0 0 0 3)) ///
replace tex nonotes label nobaselevels 


################## TABLE D.3

est clear

global controls3 "njobs curr_commute curr_startyear curr_russellgroup curr_raescore_2008 curr_acfem_pertotal2015 curr_av_fem_sal13"
global controls3_2 "mothersage_threshold nkids_threshold njobs curr_commute curr_startyear curr_russellgroup curr_raescore_2008 curr_acfem_pertotal2015 curr_av_fem_sal13"

capture noisily eststo: reghdfe curr_professor $controls3 treat_alt , vce(robust) absorb(birth_decade curr_responsibmiss curr_dept)
		estadd loc controls "Yes", replace							
		estadd loc HEI "Yes", replace
		estadd loc dept "Yes", replace							
		estadd loc age2012 "No", replace							
		
capture noisily eststo: reghdfe curr_professor $controls3_2 treat_3rd_alt , vce(robust) absorb(birth_decade curr_responsibmiss curr_dept)
		estadd loc controls "Yes", replace							
		estadd loc HEI "Yes", replace
		estadd loc dept "Yes", replace							
		estadd loc dept "Yes", replace							
		estadd loc age2012 "Yes", replace							

capture noisily eststo: reghdfe curr_gross_av_dummy $controls3 treat_alt , vce(robust) absorb(birth_decade curr_responsibmiss curr_dept)
		estadd loc controls "Yes", replace							
		estadd loc HEI "Yes", replace
		estadd loc dept "Yes", replace							
		estadd loc age2012 "No", replace							
		
capture noisily eststo: reghdfe curr_gross_av_dummy $controls3_2 treat_3rd_alt , vce(robust) absorb(birth_decade curr_responsibmiss curr_dept)
		estadd loc controls "Yes", replace							
		estadd loc HEI "Yes", replace
		estadd loc dept "Yes", replace			
		estadd loc age2012 "Yes", replace							
		
esttab est* using "Tables\Table_D3", star(* .1 ** .05 *** .01) se b(3) r2 ///
keep(treat_alt treat_3rd_alt _cons) ///
order(treat_alt treat_3rd_alt _cons) ///
varl(treat_alt "Becoming Mother (Vs. Not)" treat_3rd_alt "First Child Post-2012 (Vs. Latest Pre-2010)" _cons "Constant") ///
s(controls HEI dept age2012 N r2, ///
label("Indiv. Controls" "HEI Controls" "Research Area FEs" "Age and N. Children by 2012" "Observations" "R-squared") ///
fmt(0 0 0 0 0 3)) ///
replace tex nonotes label nobaselevels 


################## TABLE D.4

/*REDACTED FOR ANONYMITY

est clear

glo controls3 "njobs curr_commute curr_startyear curr_russellgroup curr_raescore_2008 curr_acfem_pertotal2015 curr_av_fem_sal13"
glo controls3_2 "mothersage_threshold_11 nkids_threshold_11 njobs curr_commute curr_startyear curr_russellgroup curr_raescore_2008 curr_acfem_pertotal2015 curr_av_fem_sal13"

capture noisily eststo: reghdfe curr_professor $controls3 treat_alt_11 , vce(robust) absorb(birth_decade curr_responsibmiss curr_dept)
		estadd loc controls "Yes", replace							
		estadd loc HEI "Yes", replace
		estadd loc dept "Yes", replace							
		estadd loc age2012 "No", replace							
		
capture noisily eststo: reghdfe curr_professor $controls3_2 treat_2nd_alt_11 , vce(robust) absorb(birth_decade curr_responsibmiss curr_dept)
		estadd loc controls "Yes", replace							
		estadd loc HEI "Yes", replace
		estadd loc dept "Yes", replace							
		estadd loc dept "Yes", replace							
		estadd loc age2012 "Yes", replace							

capture noisily eststo: reghdfe curr_gross_av_dummy $controls3 treat_alt_11 , vce(robust) absorb(birth_decade curr_responsibmiss curr_dept)
		estadd loc controls "Yes", replace							
		estadd loc HEI "Yes", replace
		estadd loc dept "Yes", replace							
		estadd loc age2012 "No", replace							
		
capture noisily eststo: reghdfe curr_gross_av_dummy $controls3_2 treat_2nd_alt_11 , vce(robust) absorb(birth_decade curr_responsibmiss curr_dept)
		estadd loc controls "Yes", replace							
		estadd loc HEI "Yes", replace
		estadd loc dept "Yes", replace			
		estadd loc age2012 "Yes", replace
		
esttab est* using "Tables\Table_D4", star(* .1 ** .05 *** .01) se b(3) r2 ///
keep(treat_alt_11 treat_2nd_alt_11 _cons) ///
order(treat_alt_11 treat_2nd_alt_11 _cons) ///
varl(treat_alt_11 "Becoming Mother (Vs. Not)" treat_2nd_alt_11 "Youngest Child Post-2011 (Vs. Pre-2009)" _cons "Constant") ///
s(controls HEI dept age2012 N r2, ///
label("Indiv. Controls" "HEI Controls" "Research Area FEs" "Age and N. Children by 2011" "Observations" "R-squared" "No. of HEIs") ///
fmt(0 0 0 0 0 3)) ///
replace tex nonotes label nobaselevels 
*/

################## TABLE D.5

/*REDACTED FOR ANONYMITY

est clear

glo controls3 "njobs curr_commute curr_startyear curr_russellgroup curr_raescore_2008 curr_acfem_pertotal2015 curr_av_fem_sal13"
glo controls3_2 "mothersage_threshold_13 nkids_threshold_13 njobs curr_commute curr_startyear curr_russellgroup curr_raescore_2008 curr_acfem_pertotal2015 curr_av_fem_sal13"

capture noisily eststo: reghdfe curr_professor $controls3 treat_alt_13 , vce(robust) absorb(birth_decade curr_responsibmiss curr_dept)
		estadd loc controls "Yes", replace							
		estadd loc HEI "Yes", replace
		estadd loc dept "Yes", replace							
		estadd loc age2012 "No", replace							
		
capture noisily eststo: reghdfe curr_professor $controls3_2 treat_2nd_alt_13 , vce(robust) absorb(birth_decade curr_responsibmiss curr_dept)
		estadd loc controls "Yes", replace							
		estadd loc HEI "Yes", replace
		estadd loc dept "Yes", replace							
		estadd loc dept "Yes", replace							
		estadd loc age2012 "Yes", replace							

capture noisily eststo: reghdfe curr_gross_av_dummy $controls3 treat_alt_13 , vce(robust) absorb(birth_decade curr_responsibmiss curr_dept)
		estadd loc controls "Yes", replace							
		estadd loc HEI "Yes", replace
		estadd loc dept "Yes", replace							
		estadd loc age2012 "No", replace							
		
capture noisily eststo: reghdfe curr_gross_av_dummy $controls3_2 treat_2nd_alt_13 , vce(robust) absorb(birth_decade curr_responsibmiss curr_dept)
		estadd loc controls "Yes", replace							
		estadd loc HEI "Yes", replace
		estadd loc dept "Yes", replace			
		estadd loc age2012 "Yes", replace
		
esttab est* using "Tables\Table_D5", star(* .1 ** .05 *** .01) se b(3) r2 ///
keep(treat_alt_13 treat_2nd_alt_13 _cons) ///
order(treat_alt_13 treat_2nd_alt_13 _cons) ///
varl(treat_alt_13 "Becoming Mother (Vs. Not)" treat_2nd_alt_13 "Youngest Child Post-2013 (Vs. Pre-2011)" _cons "Constant") ///
s(controls HEI dept age2012 N r2, ///
label("Indiv. Controls" "HEI Controls" "Research Area FEs" "Age and N. Children by 2013" "Observations" "R-squared" "No. of HEIs") ///
fmt(0 0 0 0 0 3)) ///
replace tex nonotes label nobaselevels 

*/

################## TABLE D.6

/*REDACTED FOR ANONYMITY

est clear

glo controls3 "njobs curr_commute curr_startyear curr_russellgroup curr_raescore_2008 curr_acfem_pertotal2015 curr_av_fem_sal13"
glo controls3_2 "mothersage_threshold_14 nkids_threshold_14 njobs curr_commute curr_startyear curr_russellgroup curr_raescore_2008 curr_acfem_pertotal2015 curr_av_fem_sal13"

capture noisily eststo: reghdfe curr_professor $controls3 treat_alt_14 , vce(robust) absorb(birth_decade curr_responsibmiss curr_dept)
		estadd loc controls "Yes", replace							
		estadd loc HEI "Yes", replace
		estadd loc dept "Yes", replace							
		estadd loc age2012 "No", replace							
		
capture noisily eststo: reghdfe curr_professor $controls3_2 treat_2nd_alt_14 , vce(robust) absorb(birth_decade curr_responsibmiss curr_dept)
		estadd loc controls "Yes", replace							
		estadd loc HEI "Yes", replace
		estadd loc dept "Yes", replace							
		estadd loc dept "Yes", replace							
		estadd loc age2012 "Yes", replace							

capture noisily eststo: reghdfe curr_gross_av_dummy $controls3 treat_alt_14 , vce(robust) absorb(birth_decade curr_responsibmiss curr_dept)
		estadd loc controls "Yes", replace							
		estadd loc HEI "Yes", replace
		estadd loc dept "Yes", replace							
		estadd loc age2012 "No", replace							
		
capture noisily eststo: reghdfe curr_gross_av_dummy $controls3_2 treat_2nd_alt_14 , vce(robust) absorb(birth_decade curr_responsibmiss curr_dept)
		estadd loc controls "Yes", replace							
		estadd loc HEI "Yes", replace
		estadd loc dept "Yes", replace			
		estadd loc age2012 "Yes", replace
		
esttab est* using "Tables\Table_D6", star(* .1 ** .05 *** .01) se b(3) r2 ///
keep(treat_alt_14 treat_2nd_alt_14 _cons) ///
order(treat_alt_14 treat_2nd_alt_14 _cons) ///
varl(treat_alt_14 "Becoming Mother (Vs. Not)" treat_2nd_alt_14 "Youngest Child Post-2014 (Vs. Pre-2012)" _cons "Constant") ///
s(controls HEI dept age2012 N r2, ///
label("Indiv. Controls" "HEI Controls" "Research Area FEs" "Age and N. Children by 2014" "Observations" "R-squared" "No. of HEIs") ///
fmt(0 0 0 0 0 3)) ///
replace tex nonotes label nobaselevels 

*/ 

################## TABLE D.7

/*REDACTED FOR ANONYMITY

est clear

glo controls3 "njobs curr_commute curr_startyear curr_russellgroup curr_raescore_2008 curr_acfem_pertotal2015 curr_av_fem_sal13"
glo controls3_2 "mothersage_threshold_15 nkids_threshold_15 njobs curr_commute curr_startyear curr_russellgroup curr_raescore_2008 curr_acfem_pertotal2015 curr_av_fem_sal13"

capture noisily eststo: reghdfe curr_professor $controls3 treat_alt_15 , vce(robust) absorb(birth_decade curr_responsibmiss curr_dept)
		estadd loc controls "Yes", replace							
		estadd loc HEI "Yes", replace
		estadd loc dept "Yes", replace							
		estadd loc age2012 "No", replace							
		
capture noisily eststo: reghdfe curr_professor $controls3_2 treat_2nd_alt_15 , vce(robust) absorb(birth_decade curr_responsibmiss curr_dept)
		estadd loc controls "Yes", replace							
		estadd loc HEI "Yes", replace
		estadd loc dept "Yes", replace							
		estadd loc dept "Yes", replace							
		estadd loc age2012 "Yes", replace							

capture noisily eststo: reghdfe curr_gross_av_dummy $controls3 treat_alt_15 , vce(robust) absorb(birth_decade curr_responsibmiss curr_dept)
		estadd loc controls "Yes", replace							
		estadd loc HEI "Yes", replace
		estadd loc dept "Yes", replace							
		estadd loc age2012 "No", replace							
		
capture noisily eststo: reghdfe curr_gross_av_dummy $controls3_2 treat_2nd_alt_15 , vce(robust) absorb(birth_decade curr_responsibmiss curr_dept)
		estadd loc controls "Yes", replace							
		estadd loc HEI "Yes", replace
		estadd loc dept "Yes", replace			
		estadd loc age2012 "Yes", replace
		
esttab est* using "Tables\Table_D7", star(* .1 ** .05 *** .01) se b(3) r2 ///
keep(treat_alt_15 treat_2nd_alt_15 _cons) ///
order(treat_alt_15 treat_2nd_alt_15 _cons) ///
varl(treat_alt_15 "Becoming Mother (Vs. Not)" treat_2nd_alt_15 "Youngest Child Post-2015 (Vs. Pre-2013)" _cons "Constant") ///
s(controls HEI dept age2012 N r2, ///
label("Indiv. Controls" "HEI Controls" "Research Area FEs" "Age and N. Children by 2015" "Observations" "R-squared" "No. of HEIs") ///
fmt(0 0 0 0 0 3)) ///
replace tex nonotes label nobaselevels 
     
*/
	 
	 

	 
	 
	 
	 
	 
	 
	 
**************************************************************************************************************************************************************************************************************************************************** CORRELATES OF PENALTY ************************************************************* 
********************************************************************************************************************************************************************************
	 
################## ################## ##################  
################## MAIN
################## ################## ################## 

################## TABLE 7

est clear

glo controls "young_mageatbirth young_permanent young_rolehigh young_rolelow young_responsiblow young_noresponsibrole"

glo uni2 "young_russellgroup young_raescore_2008 young_av_fem_sal13 young_acfem_pertotal2015" 


capture noisily eststo: reghdfe curr_gross young_p1_fwe_ml_d $controls $uni2 , vce(robust) absorb(curr_responsibmiss young_newbirthdecade birth_decade curr_dept nkids)
		estadd loc controls "Yes", replace							
		estadd loc HEI "Yes", replace
		estadd loc dept "Yes", replace			
		estadd loc decade "Yes", replace			
		estadd loc birth "Yes", replace			
		estadd loc children "Yes", replace
		
capture noisily eststo: reghdfe curr_gross young_childcare_dur_dummy young_supportive_before_table young_excluded_dummy young_relief young_lesshelp young_morehelp2 young_nohelp2 $controls $uni2 , vce(robust) absorb(curr_responsibmiss young_newbirthdecade birth_decade curr_dept nkids)
		estadd loc controls "Yes", replace							
		estadd loc HEI "Yes", replace
		estadd loc dept "Yes", replace			
		estadd loc decade "Yes", replace			
		estadd loc birth "Yes", replace			
		estadd loc children "Yes", replace
		
capture noisily eststo: reghdfe curr_professor young_p1_fwe_ml_d $controls $uni2 , vce(robust) absorb(curr_responsibmiss young_newbirthdecade birth_decade curr_dept nkids)
		estadd loc controls "Yes", replace							
		estadd loc HEI "Yes", replace
		estadd loc dept "Yes", replace			
		estadd loc decade "Yes", replace			
		estadd loc birth "Yes", replace			
		estadd loc children "Yes", replace
		
capture noisily eststo: reghdfe curr_professor young_childcare_dur_dummy young_supportive_before_table young_excluded_dummy young_relief young_lesshelp young_morehelp2 young_nohelp2 $controls $uni2 , vce(robust) absorb(curr_responsibmiss young_newbirthdecade birth_decade curr_dept nkids)
		estadd loc controls "Yes", replace							
		estadd loc HEI "Yes", replace
		estadd loc dept "Yes", replace			
		estadd loc decade "Yes", replace			
		estadd loc birth "Yes", replace			
		estadd loc children "Yes", replace
		
esttab est* using "Tables\Table_7", star(* .1 ** .05 *** .01) se b(3) r2 ///
keep(young_p1_fwe_ml_d young_childcare_dur_dummy young_supportive_before_table young_excluded_dummy young_relief young_lesshelp young_morehelp2 young_nohelp2 _cons) ///
order(young_p1_fwe_ml_d young_childcare_dur_dummy young_supportive_before_table young_excluded_dummy young_relief young_lesshelp young_morehelp2 young_nohelp2  _cons) ///
varl(young_p1_fwe_ml_d "FWE at Latest Childbirth" _cons "Constant") ///
s(controls HEI dept decade birth children N r2, ///
label("Controls at Childbirth" "HEI Features at Childbirth" "Research Area FEs" "Birth Decade FEs" "Childbirth Decade FEs" "N. Children FEs" "Observations" "R-squared") ///
fmt(0 0 0 0 0 0 0 3)) ///
replace tex nonotes label nobaselevels 


################## TABLE 8

est clear

glo controls "young_permanent young_rolehigh young_rolelow young_responsiblow young_noresponsibrole"

glo uni2 "young_russellgroup young_raescore_2008 young_av_fem_sal13 young_acfem_pertotal2015" 

capture noisily eststo: reghdfe curr_satis_reg young_p1_fwe_ml_d $controls $uni2 , vce(robust) absorb(curr_responsibmiss young_newbirthdecade birth_decade curr_dept nkids)
		estadd loc controls "Yes", replace							
		estadd loc HEI "Yes", replace
		estadd loc dept "Yes", replace			
		estadd loc decade "Yes", replace			
		estadd loc birth "Yes", replace			
		estadd loc children "Yes", replace
		
capture noisily eststo: reghdfe curr_satis_reg young_childcare_dur_dummy young_supportive_before_table young_excluded_dummy young_relief young_lesshelp young_morehelp2 young_nohelp2 $controls $uni2 , vce(robust) absorb(curr_responsibmiss young_newbirthdecade birth_decade curr_dept nkids)
		estadd loc controls "Yes", replace							
		estadd loc HEI "Yes", replace
		estadd loc dept "Yes", replace			
		estadd loc decade "Yes", replace			
		estadd loc birth "Yes", replace			
		estadd loc children "Yes", replace
		
capture noisily eststo: reghdfe curr_salfair_reg young_p1_fwe_ml_d $controls $uni2 , vce(robust) absorb(curr_responsibmiss young_newbirthdecade birth_decade curr_dept nkids)
		estadd loc controls "Yes", replace							
		estadd loc HEI "Yes", replace
		estadd loc dept "Yes", replace			
		estadd loc decade "Yes", replace			
		estadd loc birth "Yes", replace			
		estadd loc children "Yes", replace
		
capture noisily eststo: reghdfe curr_salfair_reg young_childcare_dur_dummy young_supportive_before_table young_excluded_dummy young_relief young_lesshelp young_morehelp2 young_nohelp2 $controls $uni2 , vce(robust) absorb(curr_responsibmiss young_newbirthdecade birth_decade curr_dept nkids)
		estadd loc controls "Yes", replace							
		estadd loc HEI "Yes", replace
		estadd loc dept "Yes", replace			
		estadd loc decade "Yes", replace			
		estadd loc birth "Yes", replace			
		estadd loc children "Yes", replace
		
capture noisily eststo: reghdfe curr_resfair_reg young_p1_fwe_ml_d $controls $uni2 , vce(robust) absorb(curr_responsibmiss young_newbirthdecade birth_decade curr_dept nkids)
		estadd loc controls "Yes", replace							
		estadd loc HEI "Yes", replace
		estadd loc dept "Yes", replace			
		estadd loc decade "Yes", replace			
		estadd loc birth "Yes", replace			
		estadd loc children "Yes", replace
		
capture noisily eststo: reghdfe curr_resfair_reg young_childcare_dur_dummy young_supportive_before_table young_excluded_dummy young_relief young_lesshelp young_morehelp2 young_nohelp2 $controls $uni2 , vce(robust) absorb(curr_responsibmiss young_newbirthdecade birth_decade curr_dept nkids)
		estadd loc controls "Yes", replace							
		estadd loc HEI "Yes", replace
		estadd loc dept "Yes", replace			
		estadd loc decade "Yes", replace			
		estadd loc birth "Yes", replace			
		estadd loc children "Yes", replace
		
esttab est* using "Tables\Table_8", star(* .1 ** .05 *** .01) se b(3) r2 ///
keep(young_p1_fwe_ml_d young_childcare_dur_dummy young_supportive_before_table young_excluded_dummy young_relief young_lesshelp young_morehelp2 young_nohelp2 _cons) ///
order(young_p1_fwe_ml_d young_childcare_dur_dummy young_supportive_before_table young_excluded_dummy young_relief young_lesshelp young_morehelp2 young_nohelp2  _cons) ///
varl(young_p1_fwe_ml_d "FWE at Latest Childbirth" _cons "Constant") ///
s(controls HEI dept decade birth children N r2, ///
label("Controls at Childbirth" "HEI Features at Childbirth" "Research Area FEs" "Birth Decade FEs" "Childbirth Decade FEs" "N. Children FEs" "Observations" "R-squared") ///
fmt(0 0 0 0 0 0 0 3)) ///
replace tex nonotes label nobaselevels 



################## TABLE 9

est clear

glo controls "young_permanent young_rolehigh young_rolelow young_responsiblow young_noresponsibrole"

glo uni2 "young_russellgroup young_raescore_2008 young_av_fem_sal13 young_acfem_pertotal2015" 


capture noisily eststo: reghdfe young_research_maternity young_p1_fwe_ml_d $controls $uni2 , vce(robust) absorb(curr_responsibmiss young_newbirthdecade birth_decade curr_dept nkids)
		estadd loc controls "Yes", replace							
		estadd loc HEI "Yes", replace
		estadd loc dept "Yes", replace			
		estadd loc decade "Yes", replace			
		estadd loc birth "Yes", replace			
		estadd loc children "Yes", replace			

capture noisily eststo: reghdfe young_teach_maternity young_p1_fwe_ml_d $controls $uni2 , vce(robust) absorb(curr_responsibmiss young_newbirthdecade birth_decade curr_dept nkids)
		estadd loc controls "Yes", replace							
		estadd loc HEI "Yes", replace
		estadd loc dept "Yes", replace			
		estadd loc decade "Yes", replace			
		estadd loc birth "Yes", replace			
		estadd loc children "Yes", replace			

capture noisily eststo: reghdfe young_admin_maternity young_p1_fwe_ml_d $controls $uni2 , vce(robust) absorb(curr_responsibmiss young_newbirthdecade birth_decade curr_dept nkids)
		estadd loc controls "Yes", replace							
		estadd loc HEI "Yes", replace
		estadd loc dept "Yes", replace			
		estadd loc decade "Yes", replace			
		estadd loc birth "Yes", replace			
		estadd loc children "Yes", replace			

capture noisily eststo: reghdfe young_mentor_maternity young_p1_fwe_ml_d $controls $uni2 , vce(robust) absorb(curr_responsibmiss young_newbirthdecade birth_decade curr_dept nkids)
		estadd loc controls "Yes", replace							
		estadd loc HEI "Yes", replace
		estadd loc dept "Yes", replace			
		estadd loc decade "Yes", replace			
		estadd loc birth "Yes", replace			
		estadd loc children "Yes", replace			

capture noisily eststo: reghdfe young_journal_maternity young_p1_fwe_ml_d $controls $uni2 , vce(robust) absorb(curr_responsibmiss young_newbirthdecade birth_decade curr_dept nkids)
		estadd loc controls "Yes", replace							
		estadd loc HEI "Yes", replace
		estadd loc dept "Yes", replace			
		estadd loc decade "Yes", replace			
		estadd loc birth "Yes", replace			
		estadd loc children "Yes", replace			

esttab est* using "Tables\Table_9", star(* .1 ** .05 *** .01) se b(3) r2 ///
keep(young_p1_fwe_ml_d _cons) ///
order(young_p1_fwe_ml_d _cons) ///
varl(young_p1_fwe_ml_d "FWE at Latest Childbirth" _cons "Constant") ///
s(controls HEI dept decade birth children N r2, ///
label("Controls at Childbirth" "HEI Features at Childbirth" "Research Area FEs" "Birth Decade FEs" "Childbirth Decade FEs" "N. Children FEs" "Observations" "R-squared") ///
fmt(0 0 0 0 0 0 0 3)) ///
replace tex nonotes label nobaselevels 


################## ################## ################## 
################## APPENDIX E
################## ################## ################## 

################## TABLE E.1

est clear

glo controls "young_mageatbirth young_permanent young_rolehigh young_rolelow young_responsiblow young_noresponsibrole"

glo uni2 "young_russellgroup young_raescore_2008 young_av_fem_sal13 young_acfem_pertotal2015" 

capture noisily eststo: reghdfe curr_gross young_p1_weeksfullsrep_ml_d $controls $uni2 , vce(robust) absorb(curr_responsibmiss young_newbirthdecade birth_decade curr_dept nkids)
		estadd loc controls "Yes", replace							
		estadd loc HEI "Yes", replace
		estadd loc dept "Yes", replace			
		estadd loc decade "Yes", replace			
		estadd loc birth "Yes", replace			
		estadd loc children "Yes", replace

capture noisily eststo: reghdfe curr_professor young_p1_weeksfullsrep_ml_d $controls $uni2 , vce(robust) absorb(curr_responsibmiss young_newbirthdecade birth_decade curr_dept nkids)
		estadd loc controls "Yes", replace							
		estadd loc HEI "Yes", replace
		estadd loc dept "Yes", replace			
		estadd loc decade "Yes", replace			
		estadd loc birth "Yes", replace			
		estadd loc children "Yes", replace
		
capture noisily eststo: reghdfe curr_satis_reg young_p1_weeksfullsrep_ml_d $controls $uni2 , vce(robust) absorb(curr_responsibmiss young_newbirthdecade birth_decade curr_dept nkids)
		estadd loc controls "Yes", replace							
		estadd loc HEI "Yes", replace
		estadd loc dept "Yes", replace			
		estadd loc decade "Yes", replace			
		estadd loc birth "Yes", replace			
		estadd loc children "Yes", replace
		
capture noisily eststo: reghdfe curr_salfair_reg young_p1_weeksfullsrep_ml_d $controls $uni2 , vce(robust) absorb(curr_responsibmiss young_newbirthdecade birth_decade curr_dept nkids)
		estadd loc controls "Yes", replace							
		estadd loc HEI "Yes", replace
		estadd loc dept "Yes", replace			
		estadd loc decade "Yes", replace			
		estadd loc birth "Yes", replace			
		estadd loc children "Yes", replace
		
capture noisily eststo: reghdfe curr_resfair_reg young_p1_weeksfullsrep_ml_d $controls $uni2 , vce(robust) absorb(curr_responsibmiss young_newbirthdecade birth_decade curr_dept nkids)
		estadd loc controls "Yes", replace							
		estadd loc HEI "Yes", replace
		estadd loc dept "Yes", replace			
		estadd loc decade "Yes", replace			
		estadd loc birth "Yes", replace			
		estadd loc children "Yes", replace
		
esttab est* using "Tables\Table_E1", star(* .1 ** .05 *** .01) se b(3) r2 ///
keep(young_p1_weeksfullsrep_ml_d _cons) ///
order(young_p1_weeksfullsrep_ml_d young_childcare_dur_dummy young_supportive_before_table young_excluded_dummy young_relief young_lesshelp young_morehelp2 young_nohelp2  _cons) ///
varl(young_p1_weeksfullsrep_ml_d "Full-Salary Weeks at Latest Childbirth" _cons "Constant") ///
s(controls HEI dept decade birth children N r2, ///
label("Controls at Childbirth" "HEI Features at Childbirth" "Research Area FEs" "Birth Decade FEs" "Childbirth Decade FEs" "N. Children FEs" "Observations" "R-squared") ///
fmt(0 0 0 0 0 0 0 3)) ///
replace tex label nobaselevels 


################## TABLE E.2

est clear

glo controls "young_permanent young_rolehigh young_rolelow young_responsiblow young_noresponsibrole"

glo uni2 "young_russellgroup young_raescore_2008 young_av_fem_sal13 young_acfem_pertotal2015" 


capture noisily eststo: reghdfe young_research_maternity young_p1_weeksfullsrep_ml_d $controls $uni2 , vce(robust) absorb(curr_responsibmiss young_newbirthdecade birth_decade curr_dept nkids)
		estadd loc controls "Yes", replace							
		estadd loc HEI "Yes", replace
		estadd loc dept "Yes", replace			
		estadd loc decade "Yes", replace			
		estadd loc birth "Yes", replace			
		estadd loc children "Yes", replace			

capture noisily eststo: reghdfe young_teach_maternity young_p1_weeksfullsrep_ml_d $controls $uni2 , vce(robust) absorb(curr_responsibmiss young_newbirthdecade birth_decade curr_dept nkids)
		estadd loc controls "Yes", replace							
		estadd loc HEI "Yes", replace
		estadd loc dept "Yes", replace			
		estadd loc decade "Yes", replace			
		estadd loc birth "Yes", replace			
		estadd loc children "Yes", replace			

capture noisily eststo: reghdfe young_admin_maternity young_p1_weeksfullsrep_ml_d $controls $uni2 , vce(robust) absorb(curr_responsibmiss young_newbirthdecade birth_decade curr_dept nkids)
		estadd loc controls "Yes", replace							
		estadd loc HEI "Yes", replace
		estadd loc dept "Yes", replace			
		estadd loc decade "Yes", replace			
		estadd loc birth "Yes", replace			
		estadd loc children "Yes", replace			

capture noisily eststo: reghdfe young_mentor_maternity young_p1_weeksfullsrep_ml_d $controls $uni2 , vce(robust) absorb(curr_responsibmiss young_newbirthdecade birth_decade curr_dept nkids)
		estadd loc controls "Yes", replace							
		estadd loc HEI "Yes", replace
		estadd loc dept "Yes", replace			
		estadd loc decade "Yes", replace			
		estadd loc birth "Yes", replace			
		estadd loc children "Yes", replace			

capture noisily eststo: reghdfe young_journal_maternity young_p1_weeksfullsrep_ml_d $controls $uni2 , vce(robust) absorb(curr_responsibmiss young_newbirthdecade birth_decade curr_dept nkids)
		estadd loc controls "Yes", replace							
		estadd loc HEI "Yes", replace
		estadd loc dept "Yes", replace			
		estadd loc decade "Yes", replace			
		estadd loc birth "Yes", replace			
		estadd loc children "Yes", replace			

esttab est* using "Tables\Table_E2", star(* .1 ** .05 *** .01) se b(3) r2 ///
keep(young_p1_weeksfullsrep_ml_d _cons) ///
order(young_p1_weeksfullsrep_ml_d _cons) ///
varl(young_p1_weeksfullsrep_ml_d "Full-Salary Weeks at Latest Childbirth" _cons "Constant") ///
s(controls HEI dept decade birth children N r2, ///
label("Controls at Childbirth" "HEI Features at Childbirth" "Research Area FEs" "Birth Decade FEs" "Childbirth Decade FEs" "N. Children FEs" "Observations" "R-squared") ///
fmt(0 0 0 0 0 0 0 3)) ///
replace tex nonotes label nobaselevels 


################## TABLE E.3

est clear

glo controls "young_mageatbirth young_permanent young_rolehigh young_rolelow young_responsiblow young_noresponsibrole"

glo uni2 "young_russellgroup young_raescore_2008 young_av_fem_sal13 young_acfem_pertotal2015" 


capture noisily eststo: logit curr_gross_av_dummy young_p1_fwe_ml_d $controls $uni2 i.curr_responsibmiss i.young_newbirthdecade i.birth_decade i.curr_dept i.nkids , vce(robust) 
		estadd loc controls "Yes", replace							
		estadd loc HEI "Yes", replace
		estadd loc dept "Yes", replace			
		estadd loc decade "Yes", replace			
		estadd loc birth "Yes", replace			
		estadd loc children "Yes", replace
		
capture noisily eststo: logit curr_gross_av_dummy young_childcare_dur_dummy young_supportive_before_table young_excluded_dummy young_relief young_lesshelp young_morehelp2 young_nohelp2 $controls $uni2 i.curr_responsibmiss i.young_newbirthdecade i.birth_decade i.curr_dept i.nkids , vce(robust) 
		estadd loc controls "Yes", replace							
		estadd loc HEI "Yes", replace
		estadd loc dept "Yes", replace			
		estadd loc decade "Yes", replace			
		estadd loc birth "Yes", replace			
		estadd loc children "Yes", replace
		
capture noisily eststo: logit curr_professor young_p1_fwe_ml_d $controls $uni2 i.curr_responsibmiss i.young_newbirthdecade i.birth_decade i.curr_dept i.nkids , vce(robust) 
		estadd loc controls "Yes", replace							
		estadd loc HEI "Yes", replace
		estadd loc dept "Yes", replace			
		estadd loc decade "Yes", replace			
		estadd loc birth "Yes", replace			
		estadd loc children "Yes", replace
		
capture noisily eststo: logit curr_professor young_childcare_dur_dummy young_supportive_before_table young_excluded_dummy young_relief young_lesshelp young_morehelp2 young_nohelp2 $controls $uni2 i.curr_responsibmiss i.young_newbirthdecade i.birth_decade i.curr_dept i.nkids , vce(robust) 
		estadd loc controls "Yes", replace							
		estadd loc HEI "Yes", replace
		estadd loc dept "Yes", replace			
		estadd loc decade "Yes", replace			
		estadd loc birth "Yes", replace			
		estadd loc children "Yes", replace
		
esttab est* using "Tables\Table_E3", star(* .1 ** .05 *** .01) se b(3) r2 ///
keep(young_p1_fwe_ml_d young_childcare_dur_dummy young_supportive_before_table young_excluded_dummy young_relief young_lesshelp young_morehelp2 young_nohelp2 _cons) ///
order(young_p1_fwe_ml_d young_childcare_dur_dummy young_supportive_before_table young_excluded_dummy young_relief young_lesshelp young_morehelp2 young_nohelp2  _cons) ///
varl(young_p1_fwe_ml_d "FWE at Latest Childbirth" _cons "Constant") ///
s(controls HEI dept decade birth children N , ///
label("Controls at Childbirth" "HEI Features at Childbirth" "Research Area FEs" "Birth Decade FEs" "Childbirth Decade FEs" "N. Children FEs" "Observations") ///
fmt(0 0 0 0 0 0 0 0)) ///
replace tex label nobaselevels 


################## TABLE E.4

est clear

glo controls "young_permanent young_rolehigh young_rolelow young_responsiblow young_noresponsibrole"

glo uni2 "young_russellgroup young_raescore_2008 young_av_fem_sal13 young_acfem_pertotal2015" 

capture noisily eststo: logit curr_satis_reg young_p1_fwe_ml_d $controls $uni2 i.curr_responsibmiss i.young_newbirthdecade i.birth_decade i.curr_dept i.nkids , vce(robust)
		estadd loc controls "Yes", replace							
		estadd loc HEI "Yes", replace
		estadd loc dept "Yes", replace			
		estadd loc decade "Yes", replace			
		estadd loc birth "Yes", replace			
		estadd loc children "Yes", replace
		
capture noisily eststo: logit curr_satis_reg young_childcare_dur_dummy young_supportive_before_table young_excluded_dummy young_relief young_lesshelp young_morehelp2 young_nohelp2 $controls $uni2 i.curr_responsibmiss i.young_newbirthdecade i.birth_decade i.curr_dept i.nkids , vce(robust)
		estadd loc controls "Yes", replace							
		estadd loc HEI "Yes", replace
		estadd loc dept "Yes", replace			
		estadd loc decade "Yes", replace			
		estadd loc birth "Yes", replace			
		estadd loc children "Yes", replace
		
capture noisily eststo: logit curr_salfair_reg young_p1_fwe_ml_d $controls $uni2 i.curr_responsibmiss i.young_newbirthdecade i.birth_decade i.curr_dept i.nkids , vce(robust)
		estadd loc controls "Yes", replace							
		estadd loc HEI "Yes", replace
		estadd loc dept "Yes", replace			
		estadd loc decade "Yes", replace			
		estadd loc birth "Yes", replace			
		estadd loc children "Yes", replace
		
capture noisily eststo: logit curr_salfair_reg young_childcare_dur_dummy young_supportive_before_table young_excluded_dummy young_relief young_lesshelp young_morehelp2 young_nohelp2 $controls $uni2 i.curr_responsibmiss i.young_newbirthdecade i.birth_decade i.curr_dept i.nkids , vce(robust)
		estadd loc controls "Yes", replace							
		estadd loc HEI "Yes", replace
		estadd loc dept "Yes", replace			
		estadd loc decade "Yes", replace			
		estadd loc birth "Yes", replace			
		estadd loc children "Yes", replace
		
capture noisily eststo: logit curr_resfair_reg young_p1_fwe_ml_d $controls $uni2 i.curr_responsibmiss i.young_newbirthdecade i.birth_decade i.curr_dept i.nkids , vce(robust)
		estadd loc controls "Yes", replace							
		estadd loc HEI "Yes", replace
		estadd loc dept "Yes", replace			
		estadd loc decade "Yes", replace			
		estadd loc birth "Yes", replace			
		estadd loc children "Yes", replace
		
capture noisily eststo: logit curr_resfair_reg young_childcare_dur_dummy young_supportive_before_table young_excluded_dummy young_relief young_lesshelp young_morehelp2 young_nohelp2 $controls $uni2 i.curr_responsibmiss i.young_newbirthdecade i.birth_decade i.curr_dept i.nkids , vce(robust)
		estadd loc controls "Yes", replace							
		estadd loc HEI "Yes", replace
		estadd loc dept "Yes", replace			
		estadd loc decade "Yes", replace			
		estadd loc birth "Yes", replace			
		estadd loc children "Yes", replace
		
esttab est* using "Tables\Table_E4", star(* .1 ** .05 *** .01) se b(3) r2 ///
keep(young_p1_fwe_ml_d young_childcare_dur_dummy young_supportive_before_table young_excluded_dummy young_relief young_lesshelp young_morehelp2 young_nohelp2 _cons) ///
order(young_p1_fwe_ml_d young_childcare_dur_dummy young_supportive_before_table young_excluded_dummy young_relief young_lesshelp young_morehelp2 young_nohelp2  _cons) ///
varl(young_p1_fwe_ml_d "FWE at Latest Childbirth" _cons "Constant") ///
s(controls HEI dept decade birth children N, ///
label("Controls at Childbirth" "HEI Features at Childbirth" "Research Area FEs" "Birth Decade FEs" "Childbirth Decade FEs" "N. Children FEs" "Observations") ///
fmt(0 0 0 0 0 0 0 0)) ///
replace tex label nobaselevels 


################## TABLE E.5

est clear

glo controls "young_permanent young_rolehigh young_rolelow young_responsiblow young_noresponsibrole"

glo uni2 "young_russellgroup young_raescore_2008 young_av_fem_sal13 young_acfem_pertotal2015" 


capture noisily eststo: logit young_research_maternity young_p1_fwe_ml_d $controls $uni2 i.curr_responsibmiss i.young_newbirthdecade i.birth_decade i.curr_dept i.nkids , vce(robust)
		estadd loc controls "Yes", replace							
		estadd loc HEI "Yes", replace
		estadd loc dept "Yes", replace			
		estadd loc decade "Yes", replace			
		estadd loc birth "Yes", replace			
		estadd loc children "Yes", replace			

capture noisily eststo: logit young_teach_maternity young_p1_fwe_ml_d $controls $uni2 i.curr_responsibmiss i.young_newbirthdecade i.birth_decade i.curr_dept i.nkids , vce(robust)
		estadd loc controls "Yes", replace							
		estadd loc HEI "Yes", replace
		estadd loc dept "Yes", replace			
		estadd loc decade "Yes", replace			
		estadd loc birth "Yes", replace			
		estadd loc children "Yes", replace			

capture noisily eststo: logit young_admin_maternity young_p1_fwe_ml_d $controls $uni2 i.curr_responsibmiss i.young_newbirthdecade i.birth_decade i.curr_dept i.nkids , vce(robust)
		estadd loc controls "Yes", replace							
		estadd loc HEI "Yes", replace
		estadd loc dept "Yes", replace			
		estadd loc decade "Yes", replace			
		estadd loc birth "Yes", replace			
		estadd loc children "Yes", replace			

capture noisily eststo: logit young_mentor_maternity young_p1_fwe_ml_d $controls $uni2 i.curr_responsibmiss i.young_newbirthdecade i.birth_decade i.curr_dept i.nkids , vce(robust)
		estadd loc controls "Yes", replace							
		estadd loc HEI "Yes", replace
		estadd loc dept "Yes", replace			
		estadd loc decade "Yes", replace			
		estadd loc birth "Yes", replace			
		estadd loc children "Yes", replace			

capture noisily eststo: logit young_journal_maternity young_p1_fwe_ml_d $controls $uni2 i.curr_responsibmiss i.young_newbirthdecade i.birth_decade i.curr_dept i.nkids , vce(robust)
		estadd loc controls "Yes", replace							
		estadd loc HEI "Yes", replace
		estadd loc dept "Yes", replace			
		estadd loc decade "Yes", replace			
		estadd loc birth "Yes", replace			
		estadd loc children "Yes", replace			

esttab est* using "Tables\Table_E5", star(* .1 ** .05 *** .01) se b(3) r2 ///
keep(young_p1_fwe_ml_d _cons) ///
order(young_p1_fwe_ml_d _cons) ///
varl(young_p1_fwe_ml_d "FWE at Latest Childbirth" _cons "Constant") ///
s(controls HEI dept decade birth children N, ///
label("Controls at Childbirth" "HEI Features at Childbirth" "Research Area FEs" "Birth Decade FEs" "Childbirth Decade FEs" "N. Children FEs" "Observations") ///
fmt(0 0 0 0 0 0 0)) ///
replace tex label nobaselevels 


################## TABLE E.6

est clear

glo controls "oldest_mageatbirth oldest_permanent oldest_rolehigh oldest_rolelow oldest_responsiblow oldest_noresponsibrole"

glo uni2 "oldest_russellgroup oldest_raescore_2008 oldest_av_fem_sal13 oldest_acfem_pertotal2015" 


capture noisily eststo: reghdfe curr_gross oldest_p1_fwe_ml_d $controls $uni2 , vce(robust) absorb(curr_responsibmiss oldest_newbirthdecade birth_decade curr_dept nkids)
		estadd loc controls "Yes", replace							
		estadd loc HEI "Yes", replace
		estadd loc dept "Yes", replace			
		estadd loc decade "Yes", replace			
		estadd loc birth "Yes", replace			
		estadd loc children "Yes", replace
		
capture noisily eststo: reghdfe curr_gross oldest_childcare_dur_dummy oldest_supportive_before_table oldest_excluded_dummy oldest_relief oldest_lesshelp oldest_morehelp2 oldest_nohelp2 $controls $uni2 , vce(robust) absorb(curr_responsibmiss oldest_newbirthdecade birth_decade curr_dept nkids)
		estadd loc controls "Yes", replace							
		estadd loc HEI "Yes", replace
		estadd loc dept "Yes", replace			
		estadd loc decade "Yes", replace			
		estadd loc birth "Yes", replace			
		estadd loc children "Yes", replace
		
capture noisily eststo: reghdfe curr_professor oldest_p1_fwe_ml_d $controls $uni2 , vce(robust) absorb(curr_responsibmiss oldest_newbirthdecade birth_decade curr_dept nkids)
		estadd loc controls "Yes", replace							
		estadd loc HEI "Yes", replace
		estadd loc dept "Yes", replace			
		estadd loc decade "Yes", replace			
		estadd loc birth "Yes", replace			
		estadd loc children "Yes", replace
		
capture noisily eststo: reghdfe curr_professor oldest_childcare_dur_dummy oldest_supportive_before_table oldest_excluded_dummy oldest_relief oldest_lesshelp oldest_morehelp2 oldest_nohelp2 $controls $uni2 , vce(robust) absorb(curr_responsibmiss oldest_newbirthdecade birth_decade curr_dept nkids)
		estadd loc controls "Yes", replace							
		estadd loc HEI "Yes", replace
		estadd loc dept "Yes", replace			
		estadd loc decade "Yes", replace			
		estadd loc birth "Yes", replace			
		estadd loc children "Yes", replace
		
esttab est* using "Tables\Table_E6", star(* .1 ** .05 *** .01) se b(3) r2 ///
keep(oldest_p1_fwe_ml_d oldest_childcare_dur_dummy oldest_supportive_before_table oldest_excluded_dummy oldest_relief oldest_lesshelp oldest_morehelp2 oldest_nohelp2 _cons) ///
order(oldest_p1_fwe_ml_d oldest_childcare_dur_dummy oldest_supportive_before_table oldest_excluded_dummy oldest_relief oldest_lesshelp oldest_morehelp2 oldest_nohelp2  _cons) ///
varl(oldest_p1_fwe_ml_d "FWE at First Childbirth" _cons "Constant") ///
s(controls HEI dept decade birth children N r2, ///
label("Controls at Childbirth" "HEI Features at Childbirth" "Research Area FEs" "Birth Decade FEs" "Childbirth Decade FEs" "N. Children FEs" "Observations" "R-squared") ///
fmt(0 0 0 0 0 0 0 3)) ///
replace tex label nobaselevels 


################## TABLE E.7

est clear

glo controls "oldest_permanent oldest_rolehigh oldest_rolelow oldest_responsiblow oldest_noresponsibrole"

glo uni2 "oldest_russellgroup oldest_raescore_2008 oldest_av_fem_sal13 oldest_acfem_pertotal2015" 

capture noisily eststo: reghdfe curr_satis_reg oldest_p1_fwe_ml_d $controls $uni2 , vce(robust) absorb(curr_responsibmiss oldest_newbirthdecade birth_decade curr_dept nkids)
		estadd loc controls "Yes", replace							
		estadd loc HEI "Yes", replace
		estadd loc dept "Yes", replace			
		estadd loc decade "Yes", replace			
		estadd loc birth "Yes", replace			
		estadd loc children "Yes", replace
		
capture noisily eststo: reghdfe curr_satis_reg oldest_childcare_dur_dummy oldest_supportive_before_table oldest_excluded_dummy oldest_relief oldest_lesshelp oldest_morehelp2 oldest_nohelp2 $controls $uni2 , vce(robust) absorb(curr_responsibmiss oldest_newbirthdecade birth_decade curr_dept nkids)
		estadd loc controls "Yes", replace							
		estadd loc HEI "Yes", replace
		estadd loc dept "Yes", replace			
		estadd loc decade "Yes", replace			
		estadd loc birth "Yes", replace			
		estadd loc children "Yes", replace
		
capture noisily eststo: reghdfe curr_salfair_reg oldest_p1_fwe_ml_d $controls $uni2 , vce(robust) absorb(curr_responsibmiss oldest_newbirthdecade birth_decade curr_dept nkids)
		estadd loc controls "Yes", replace							
		estadd loc HEI "Yes", replace
		estadd loc dept "Yes", replace			
		estadd loc decade "Yes", replace			
		estadd loc birth "Yes", replace			
		estadd loc children "Yes", replace
		
capture noisily eststo: reghdfe curr_salfair_reg oldest_childcare_dur_dummy oldest_supportive_before_table oldest_excluded_dummy oldest_relief oldest_lesshelp oldest_morehelp2 oldest_nohelp2 $controls $uni2 , vce(robust) absorb(curr_responsibmiss oldest_newbirthdecade birth_decade curr_dept nkids)
		estadd loc controls "Yes", replace							
		estadd loc HEI "Yes", replace
		estadd loc dept "Yes", replace			
		estadd loc decade "Yes", replace			
		estadd loc birth "Yes", replace			
		estadd loc children "Yes", replace
		
capture noisily eststo: reghdfe curr_resfair_reg oldest_p1_fwe_ml_d $controls $uni2 , vce(robust) absorb(curr_responsibmiss oldest_newbirthdecade birth_decade curr_dept nkids)
		estadd loc controls "Yes", replace							
		estadd loc HEI "Yes", replace
		estadd loc dept "Yes", replace			
		estadd loc decade "Yes", replace			
		estadd loc birth "Yes", replace			
		estadd loc children "Yes", replace
		
capture noisily eststo: reghdfe curr_resfair_reg oldest_childcare_dur_dummy oldest_supportive_before_table oldest_excluded_dummy oldest_relief oldest_lesshelp oldest_morehelp2 oldest_nohelp2 $controls $uni2 , vce(robust) absorb(curr_responsibmiss oldest_newbirthdecade birth_decade curr_dept nkids)
		estadd loc controls "Yes", replace							
		estadd loc HEI "Yes", replace
		estadd loc dept "Yes", replace			
		estadd loc decade "Yes", replace			
		estadd loc birth "Yes", replace			
		estadd loc children "Yes", replace
		
esttab est* using "Tables\Table_E7", star(* .1 ** .05 *** .01) se b(3) r2 ///
keep(oldest_p1_fwe_ml_d oldest_childcare_dur_dummy oldest_supportive_before_table oldest_excluded_dummy oldest_relief oldest_lesshelp oldest_morehelp2 oldest_nohelp2 _cons) ///
order(oldest_p1_fwe_ml_d oldest_childcare_dur_dummy oldest_supportive_before_table oldest_excluded_dummy oldest_relief oldest_lesshelp oldest_morehelp2 oldest_nohelp2  _cons) ///
varl(oldest_p1_fwe_ml_d "FWE at First Childbirth" _cons "Constant") ///
s(controls HEI dept decade birth children N r2, ///
label("Controls at Childbirth" "HEI Features at Childbirth" "Research Area FEs" "Birth Decade FEs" "Childbirth Decade FEs" "N. Children FEs" "Observations" "R-squared") ///
fmt(0 0 0 0 0 0 0 3)) ///
replace tex label nobaselevels 



################## TABLE E.8

est clear

glo controls "oldest_permanent oldest_rolehigh oldest_rolelow oldest_responsiblow oldest_noresponsibrole"

glo uni2 "oldest_russellgroup oldest_raescore_2008 oldest_av_fem_sal13 oldest_acfem_pertotal2015" 


capture noisily eststo: reghdfe oldest_research_maternity oldest_p1_fwe_ml_d $controls $uni2 , vce(robust) absorb(curr_responsibmiss oldest_newbirthdecade birth_decade curr_dept nkids)
		estadd loc controls "Yes", replace							
		estadd loc HEI "Yes", replace
		estadd loc dept "Yes", replace			
		estadd loc decade "Yes", replace			
		estadd loc birth "Yes", replace			
		estadd loc children "Yes", replace			

capture noisily eststo: reghdfe oldest_teach_maternity oldest_p1_fwe_ml_d $controls $uni2 , vce(robust) absorb(curr_responsibmiss oldest_newbirthdecade birth_decade curr_dept nkids)
		estadd loc controls "Yes", replace							
		estadd loc HEI "Yes", replace
		estadd loc dept "Yes", replace			
		estadd loc decade "Yes", replace			
		estadd loc birth "Yes", replace			
		estadd loc children "Yes", replace			

capture noisily eststo: reghdfe oldest_admin_maternity oldest_p1_fwe_ml_d $controls $uni2 , vce(robust) absorb(curr_responsibmiss oldest_newbirthdecade birth_decade curr_dept nkids)
		estadd loc controls "Yes", replace							
		estadd loc HEI "Yes", replace
		estadd loc dept "Yes", replace			
		estadd loc decade "Yes", replace			
		estadd loc birth "Yes", replace			
		estadd loc children "Yes", replace			

capture noisily eststo: reghdfe oldest_mentor_maternity oldest_p1_fwe_ml_d $controls $uni2 , vce(robust) absorb(curr_responsibmiss oldest_newbirthdecade birth_decade curr_dept nkids)
		estadd loc controls "Yes", replace							
		estadd loc HEI "Yes", replace
		estadd loc dept "Yes", replace			
		estadd loc decade "Yes", replace			
		estadd loc birth "Yes", replace			
		estadd loc children "Yes", replace			

capture noisily eststo: reghdfe oldest_journal_maternity oldest_p1_fwe_ml_d $controls $uni2 , vce(robust) absorb(curr_responsibmiss oldest_newbirthdecade birth_decade curr_dept nkids)
		estadd loc controls "Yes", replace							
		estadd loc HEI "Yes", replace
		estadd loc dept "Yes", replace			
		estadd loc decade "Yes", replace			
		estadd loc birth "Yes", replace			
		estadd loc children "Yes", replace			

esttab est* using "Tables\Table_E8", star(* .1 ** .05 *** .01) se b(3) r2 ///
keep(oldest_p1_fwe_ml_d _cons) ///
order(oldest_p1_fwe_ml_d _cons) ///
varl(oldest_p1_fwe_ml_d "FWE at First Childbirth" _cons "Constant") ///
s(controls HEI dept decade birth children N r2, ///
label("Controls at Childbirth" "HEI Features at Childbirth" "Research Area FEs" "Birth Decade FEs" "Childbirth Decade FEs" "N. Children FEs" "Observations" "R-squared") ///
fmt(0 0 0 0 0 0 0 3)) ///
replace tex label nobaselevels 



################## TABLE E.9

/*REDACTED FOR ANONYMITY

est clear

glo controls "oldest_mageatbirth oldest_permanent oldest_rolehigh oldest_rolelow oldest_responsiblow oldest_noresponsibrole"

glo uni2 "oldest_russellgroup oldest_raescore_2008 oldest_av_fem_sal13 oldest_acfem_pertotal2015" 


capture noisily eststo: reghdfe curr_gross fwe_avg $controls $uni2 , vce(robust) absorb(curr_responsibmiss oldest_newbirthdecade birth_decade curr_dept nkids)
		estadd loc controls "Yes", replace							
		estadd loc HEI "Yes", replace
		estadd loc dept "Yes", replace			
		estadd loc decade "Yes", replace			
		estadd loc birth "Yes", replace			
		estadd loc children "Yes", replace
		
capture noisily eststo: reghdfe curr_gross childcare_avg supportive_avg excluded_avg relief_avg lesshelp2_avg morehelp2_avg nohelp2_avg $controls $uni2 , vce(robust) absorb(curr_responsibmiss oldest_newbirthdecade birth_decade curr_dept nkids)
		estadd loc controls "Yes", replace							
		estadd loc HEI "Yes", replace
		estadd loc dept "Yes", replace			
		estadd loc decade "Yes", replace			
		estadd loc birth "Yes", replace			
		estadd loc children "Yes", replace
		
capture noisily eststo: reghdfe curr_professor fwe_avg $controls $uni2 , vce(robust) absorb(curr_responsibmiss oldest_newbirthdecade birth_decade curr_dept nkids)
		estadd loc controls "Yes", replace							
		estadd loc HEI "Yes", replace
		estadd loc dept "Yes", replace			
		estadd loc decade "Yes", replace			
		estadd loc birth "Yes", replace			
		estadd loc children "Yes", replace
		
capture noisily eststo: reghdfe curr_professor childcare_avg supportive_avg excluded_avg relief_avg lesshelp2_avg morehelp2_avg nohelp2_avg $controls $uni2 , vce(robust) absorb(curr_responsibmiss oldest_newbirthdecade birth_decade curr_dept nkids)
		estadd loc controls "Yes", replace							
		estadd loc HEI "Yes", replace
		estadd loc dept "Yes", replace			
		estadd loc decade "Yes", replace			
		estadd loc birth "Yes", replace			
		estadd loc children "Yes", replace
		
esttab est* using "Tables\Table_E9", star(* .1 ** .05 *** .01) se b(3) r2 ///
keep(fwe_avg childcare_avg supportive_avg excluded_avg relief_avg lesshelp2_avg morehelp2_avg nohelp2_avg _cons) ///
order(fwe_avg childcare_avg supportive_avg excluded_avg relief_avg lesshelp2_avg morehelp2_avg nohelp2_avg  _cons) ///
varl(fwe_avg "FWE (Avg)" _cons "Constant") ///
s(controls HEI dept decade birth children N r2, ///
label("Controls at Childbirth" "HEI Features at Childbirth" "Research Area FEs" "Birth Decade FEs" "Childbirth Decade FEs" "N. Children FEs" "Observations" "R-squared") ///
fmt(0 0 0 0 0 0 0 3)) ///
replace tex label nobaselevels 
*/


################## TABLE E.10

/*REDACTED FOR ANONYMITY

est clear

glo controls "oldest_permanent oldest_rolehigh oldest_rolelow oldest_responsiblow oldest_noresponsibrole"

glo uni2 "oldest_russellgroup oldest_raescore_2008 oldest_av_fem_sal13 oldest_acfem_pertotal2015" 

capture noisily eststo: reghdfe curr_satis_reg fwe_avg $controls $uni2 , vce(robust) absorb(curr_responsibmiss oldest_newbirthdecade birth_decade curr_dept nkids)
		estadd loc controls "Yes", replace							
		estadd loc HEI "Yes", replace
		estadd loc dept "Yes", replace			
		estadd loc decade "Yes", replace			
		estadd loc birth "Yes", replace			
		estadd loc children "Yes", replace
		
capture noisily eststo: reghdfe curr_satis_reg childcare_avg supportive_avg excluded_avg relief_avg lesshelp2_avg morehelp2_avg nohelp2_avg $controls $uni2 , vce(robust) absorb(curr_responsibmiss oldest_newbirthdecade birth_decade curr_dept nkids)
		estadd loc controls "Yes", replace							
		estadd loc HEI "Yes", replace
		estadd loc dept "Yes", replace			
		estadd loc decade "Yes", replace			
		estadd loc birth "Yes", replace			
		estadd loc children "Yes", replace
		
capture noisily eststo: reghdfe curr_salfair_reg fwe_avg $controls $uni2 , vce(robust) absorb(curr_responsibmiss oldest_newbirthdecade birth_decade curr_dept nkids)
		estadd loc controls "Yes", replace							
		estadd loc HEI "Yes", replace
		estadd loc dept "Yes", replace			
		estadd loc decade "Yes", replace			
		estadd loc birth "Yes", replace			
		estadd loc children "Yes", replace
		
capture noisily eststo: reghdfe curr_salfair_reg childcare_avg supportive_avg excluded_avg relief_avg lesshelp2_avg morehelp2_avg nohelp2_avg $controls $uni2 , vce(robust) absorb(curr_responsibmiss oldest_newbirthdecade birth_decade curr_dept nkids)
		estadd loc controls "Yes", replace							
		estadd loc HEI "Yes", replace
		estadd loc dept "Yes", replace			
		estadd loc decade "Yes", replace			
		estadd loc birth "Yes", replace			
		estadd loc children "Yes", replace
		
capture noisily eststo: reghdfe curr_resfair_reg fwe_avg $controls $uni2 , vce(robust) absorb(curr_responsibmiss oldest_newbirthdecade birth_decade curr_dept nkids)
		estadd loc controls "Yes", replace							
		estadd loc HEI "Yes", replace
		estadd loc dept "Yes", replace			
		estadd loc decade "Yes", replace			
		estadd loc birth "Yes", replace			
		estadd loc children "Yes", replace
		
capture noisily eststo: reghdfe curr_resfair_reg childcare_avg supportive_avg excluded_avg relief_avg lesshelp2_avg morehelp2_avg nohelp2_avg $controls $uni2 , vce(robust) absorb(curr_responsibmiss oldest_newbirthdecade birth_decade curr_dept nkids)
		estadd loc controls "Yes", replace							
		estadd loc HEI "Yes", replace
		estadd loc dept "Yes", replace			
		estadd loc decade "Yes", replace			
		estadd loc birth "Yes", replace			
		estadd loc children "Yes", replace
		
esttab est* using "Tables\Table_E10", star(* .1 ** .05 *** .01) se b(3) r2 ///
keep(fwe_avg childcare_avg supportive_avg excluded_avg relief_avg lesshelp2_avg morehelp2_avg nohelp2_avg _cons) ///
order(fwe_avg childcare_avg supportive_avg excluded_avg relief_avg lesshelp2_avg morehelp2_avg nohelp2_avg  _cons) ///
varl(fwe_avg "FWE (Avg)" _cons "Constant") ///
s(controls HEI dept decade birth children N r2, ///
label("Controls at Childbirth" "HEI Features at Childbirth" "Research Area FEs" "Birth Decade FEs" "Childbirth Decade FEs" "N. Children FEs" "Observations" "R-squared") ///
fmt(0 0 0 0 0 0 0 3)) ///
replace tex label nobaselevels 

*/


################## TABLE E.11

est clear

glo controls "young_permanent young_rolehigh young_rolelow young_responsiblow young_noresponsibrole"

glo uni2 "young_russellgroup young_raescore_2008 young_av_fem_sal13 young_acfem_pertotal2015" 

capture noisily eststo: reghdfe curr_satis_reg_alt young_p1_fwe_ml_d $controls $uni2 , vce(robust) absorb(curr_responsibmiss young_newbirthdecade birth_decade curr_dept nkids)
		estadd loc controls "Yes", replace							
		estadd loc HEI "Yes", replace
		estadd loc dept "Yes", replace			
		estadd loc decade "Yes", replace			
		estadd loc birth "Yes", replace			
		estadd loc children "Yes", replace
		
capture noisily eststo: reghdfe curr_satis_reg_alt young_childcare_dur_dummy young_supportive_before_table young_excluded_dummy young_relief young_lesshelp young_morehelp2 young_nohelp2 $controls $uni2 , vce(robust) absorb(curr_responsibmiss young_newbirthdecade birth_decade curr_dept nkids)
		estadd loc controls "Yes", replace							
		estadd loc HEI "Yes", replace
		estadd loc dept "Yes", replace			
		estadd loc decade "Yes", replace			
		estadd loc birth "Yes", replace			
		estadd loc children "Yes", replace
		
esttab est* using "Tables\Table_E11", star(* .1 ** .05 *** .01) se b(3) r2 ///
keep(young_p1_fwe_ml_d young_childcare_dur_dummy young_supportive_before_table young_excluded_dummy young_relief young_lesshelp young_morehelp2 young_nohelp2 _cons) ///
order(young_p1_fwe_ml_d young_childcare_dur_dummy young_supportive_before_table young_excluded_dummy young_relief young_lesshelp young_morehelp2 young_nohelp2  _cons) ///
varl(young_p1_fwe_ml_d "FWE at Latest Childbirth" _cons "Constant") ///
s(controls HEI dept decade birth children N r2, ///
label("Controls at Childbirth" "HEI Features at Childbirth" "Research Area FEs" "Birth Decade FEs" "Childbirth Decade FEs" "N. Children FEs" "Observations" "R-squared") ///
fmt(0 0 0 0 0 0 0 3)) ///
replace tex nonotes label nobaselevels 



